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Asymptotically free gauge theories of the strong interactions are constructed and analyzed. The reasons 
for doing this are recounted, including a review of renormalization-group techniques and their 
application to scaling phenomena. The renormalization-group equations are derived for Yang-Mills 
theories. The parameters that enter into the equations are calculated to lowest order and it is shown 
that these theories are asymptotically free. More specifically the effective coupling constant, which 
determines the ultraviolet behavior of the theory, vanishes for large spacelike momenta. Fermions are 
incorporated and the construction of realistic models is discussed. We propose that the strong 
interactions be mediated by a “color” gauge group which commutes with SU(3) x SU(3). The problem 
of symmetry breaking is discussed. It appears likely that this would have a dynamical origin. It 1s 
suggested that the gauge symmetry might not be broken and that the severe infrared singularities 
prevent the occurrence of noncolor singlet physical states. The deep-inelastic structure functions, as well 
as the electron-positron total annihilation cross section are analyzed. Scaling obtains up to calculable 
logarithmic corrections, and the naive light-cone or parton-model results follow. The problems of 
incorporating scalar mesons and breaking the symmetry by the Higgs mechanism are explained in detail. 


I. INTRODUCTION 


In recent years the renormalization group has 
played an increasingly important role in the study 
of the asymptotic behavior of renormalizable field 
theories. This approach has acquired new impor- 
tance due to the recent discovery that non-Abelian 
gauge theories are asymptotically free.*™ In this 
paper we shall amplify and extend the results re- 
ported in Ref. 1. 

The renormalization group dates from the funda- 
mental work of Gell-Mann and Low,* who studied 
the asymptotic behavior of the photon propagator 
in quantum electrodynamics. The remarkable 
discovery of Gell-Mann and Low was that the as- 
ymptotic form of the photon propagator was deter- 
mined by the zeros of a certain, calculable, func- 
tion of the coupling constant and not by the actual 
value of the charge. The renormalization-group 
equations were extended by Bogolubov and Shirkov 
to the vertex function® and employed to analyze the 
ultraviolet and infrared behavior of quantum elec- 
trodynamics and other field theories.” (For a 
review of this work see Ref. 8.) 

The basic idea underlying the renormalization- 
group equations is very simple. A renormalizable 
field theory contains two types of parameters— 
masses or coupling constants with positive dimen- 
sions of mass (i.e., due to Myy or àg? terms in 
the Lagrangian) and dimensionless coupling con- 
stants (i.e., due to Af or yr" PA, terms in the 
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Lagrangian). Coupling constants with negative di- 
mensions of mass give nonrenormalizable theo- 
ries. If one considers a Green’s function for 
large and spacelike momenta (so as to exclude any 
Landau singularities), then one would expect that 
the generalized mass terms in the Lagrangian 
(Myy or Ap?) could be neglected. In other words 
the leading asymptotic behavior of the Green’s 
functions should be the same as would be calculated 
ina massless theory. This can be proved, to any 
finite order in perturbation theory, by using 
Weinberg’s theorem.” The massless theory con- 
tains no dimensional parameters to set the scale 
of momenta, therefore one might expect that the 
asymptotic behavior of the amplitudes would be 
determined by pure dimensional analysis. This is 
called naive or cannonical scaling. It does not 
occur in practice, since the massless theory does 
contain a hidden dimensional parameter. This 
parameter, u, must be introduced in order to 
perform the subtractions necessary to renormalize 
the theory and render it finite. Due to infrared 
singularities these subtractions, for the massless 
theory, must be performed off shell, say at some 
spacelike momenta p° =-— u°’. The subtractions 
then define the physical coupling constants and the 
scale of the fields (which are determined by the 
wave-function renormalization constants). The 
subtraction point, u, is arbitrary. If we change 
the subtraction point the net effect is to change the 
value of the coupling constants and the scale of the 
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fields. This fact is expressed by the renormaliza- 
tion-group equations. 

A change in the subtraction point, u, is equiva - 
lent to a change in the scale of all momenta since 
u is the sole parameter that fixes the momentum 
scale. Therefore one can use the renormalization 
group to relate Green’s functions for one set of 
momenta and coupling constant to Green’s function 
for a scaled set of momenta and different values 
of the coupling constants. To do this one only 
needs to know the functions which determine the 
change in coupling constant and scale of fields due 
to a change in u. In particular the asymptotic 
form of the amplitudes (for large spacelike mo- 
menta) can be related to amplitudes at some fixed 
momenta and an effective coupling constant (which 
would be determined by performing the subtrac- 
tions at asymptotic values of the momenta). The 
asymptotic value of this effective coupling con- 
stant will determine the ultraviolet behavior of 
the theory. It is given by the zeros of a calculable 
function, which are called fixed points of the re- 
normalization group. 

For about ten years there was little interest in 
the renormalization group.’° This was probably 
due to the following reasons. First, the renormal- 
ization group provided information about Green’s 
functions for large spacelike momenta, which are 
of no direct physical interest. This approach is 
much less informative about Minkowski momenta 
and on-shell amplitudes—since we lack an exten- 
sion to this region of Weinbereg’s powerful theo- 
rem.” Second, it was soon discovered, in all the 
cases investigated at the time, that the ultraviolet 
behavior was not calculable using perturbation 
theory. Thus it appeared that the renormalization 
group provided a framework in which one could 
discuss, but not calculate, the asymptotic behavior 
of amplitudes in a physically uninteresting region. 

This situation has changed in the last few years 
due to the following developments. First there 
was an increased interest in the matrix elements 
of local currents at short distances” (or their 
Fourier transforms for large spacelike momenta). 
This interest was further increased by the advent 
of new experiments on deep-inelastic lepton-had- 
ron scattering at SLAC, and by Bjorken’s predic- 
tion and the experimental indications of scaling.” 
The theoretical framework for the discussion of 
products of currents at short distances was pro- 
vided by Wilson’s operator-product expansion.** 
Although Wilson had emphasized that, in general, 
one does not expect naive scaling, the experimental 
indications of Bjorken scaling motivated the de- 
velopment of phenomenological scaling models. 
One approach was to abstract from free-field the- 
ory, or interacting-field theories with an ultra- 


violet cutoff, the short-distance structure of cur- 
rent products, the other was to hypothesize the 
existence of pointlike constituents of hadrons 
(partons).'? Assuming a quark field -theoretic mod- 
el, or that the partons had the quantum numbers 

of quarks, many relations and sum rules were de- 
rived. At present the experimental data are con- 
sistent with spin-3 (Ref. 16) quarklike constituents. 
All attempts to provide dynamical explanations for 
these models were unsuccessful. 

Meanwhile significant developments of the re- 
normalization group techniques were made. The 
equations that determine the change in momentum 
scale of massive field theories were derived by 
Callan” and Symanzik.® These equations, for 
large spacelike momenta, reduce to the differen- 
tial form of the renormalization-group equations 
already derived by Orsiannikov in a little-noticed 
paper in 1956."° In addition it was realized that 
the renormalization-group approach is the key to 
asymptotic behavior of the coefficient functions in 
Wilson’s operator-product expansion, and thus the 
related behavior of electroproduction structure 
functions in the Bjorken region.”°~?8 

In particular it was realized that Bjorken scaling 
could be understood within the framework of the 
renormalization group if there was an ultraviolet- 
(UV) stable fixed point of the renormalization 
group. At the fixed point, however, the anomalous 
dimensions of the relevant operators in Wilson’s 
expansion would all have to vanish. All the indica- 
tions are that this can only occur if the value of the 
fixed point of the renormalization group is zero.?*°?° 
This has recently been proved in a large class of 
field theories.*® In such a theory the effective 
coupling constant vanishes for large spacelike mo- 
menta and we describe this phenomena by saying 
that the theory is asymptotically free. An asymp- 
totically free theory will exhibit Bjorken scaling 
(up to, perhaps, logarithmic corrections) and in 
addition will lead to all the naive light-cone- or 
parton-model results.”° 

The possibility that a given field theory is as- 
ymptotically free is easily explored by simple 
perturbation-theory calculations. Quantum elec- 
trodynamics was known, from the original work of 
Gell-Mann and Low, not to be asymptotically free. 
Zee extended this result to scalar-fermion theo- 
ries involving one coupling constant.?” Recently 
Coleman and one of us (D.J.G.) have proved that 
no theory which does not involve non-Abelian 
gauge mesons can be asymptotically free.” 
Together with the recent discovery that non- 
Abelian gauge theories are asymptotically freet” 
these developments provide a compelling case for 
a non-Abelian gauge theory of the strong interac- 
tions. Indeed if one accepts the renormalization- 
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group approach and the experimental reality of 
Bjorken scaling as an asymptotic phenomenon then 
there is, probably, no other choice.” This pos- 
sibility is explored in the following. 

In Sec. II we shall outline the derivation of the 
renormalization-group equations for pure Yang- 
Mills theories (which involve only gauge fields). 
These equations are discussed in greater detail in 
Appendix A. We also discuss the notion of the ef- 
fective coupling constant and exhibit the solution 
of the renormalization-group equations. 

In Sec. III we calculate the renormalization- 
group parameters for pure Yang-Mills theories. 

In Sec. IV we incorporate fermions into the gauge 
theories, without destroying their asymptotic 
freedom. The large-momentum behavior of the 
effective coupling constant and Green’s functions 
is derived. | | 

In Sec. V the construction of realistic physical 
models of the strong interactions is discussed. 
We analyze the structure functions of deep-inelas- 
tic scattering and the total electron- positron an- 
nihilation cross section in these models. These 
applications will be explored further in a forth- 
coming publication.*° The major problem remain- 
ing in these gauge theories is how to break the 
gauge symmetry and provide masses for the vector 
mesons. Various dynamical possibilities are dis- 
cussed in Sec. V. 

In Sec. VI we incorporate scalar mesons into 
asymptotically free gauge theories. The difficul- 
ties encountered in achieving this are described, 
and the failure to construct models in which the 
Higgs mechanism generates masses for all the 
vector mesons is explained. 

Section VII contains some concluding remarks. 


H. THE RENORMALIZATION-GROUP EQUATIONS 
FOR YANG-MILLS THEORIES 


In this section we shall derive the renormaliza- 
tion-group equations for a pure Yang-Mills theory. 
The only restriction on the form of the theory will 
be the requirement that the gauge group be non- 
Abelian and semisimple. 

The classical Yang-Mills Lagrangian density is 

==3Tr{(8,B,-8,B,— g[B, B,]|)*}, (2.1) 
where 
B (x) =B Ota (2.2) 


is a matrix of Hermitian vector fields (summation 
over repeated indices is implied). The matrices 
t, generate a semisimple Lie group G: 


[tastu] = iCisete (2.3) 


and are normalized according to 


TEIL laO (2.4) 
b a 


This Lagrangian is singular due to its invariance 
under the gauge group. Therefore a proper quan- 
tization of (2.1) necessitates the addition of a 
gauge-fixing term to the Lagrangian, say 
-(1/a)Tr{(ə" B ,}}. The presence of this term 
then requires the appearance of Feynman- Faddeev- 
Popov ghosts.** The net result is that the effective 
Lagrangian used to derive the Feynman rules is 


=-3 Tr{(0,B,- 9,B, - g[ By B,] Y} 


1 
q 


+2 Tr{a ġ*ð p- gə >*[ By, l}, (2.5) 


where ọ= °t, is a massless, complex, scalar 
field which propagates in closed loops only and 

obeys Fermi statistics. The resulting Feynman 
rules are summarized in Fig. 1. 

Due to the presence of massless particles and 
the resulting infrared singularities it would appear 
that the S matrix does not exist, at least in pertur- 
bation theory, for the pure Yang-Mills theory. 

One can however consider the off-shell Green’s 
function for such theories, at all but exceptional 
momenta.?? In the following, this restriction is 
always to be understood. The large-momentum 
behavior of such Green’s functions is not without 
physical interest since it bears directly on the 
ultraviolet behavior of more realistic models with 
symmetry-breaking and mass terms. 

Up till now we have discussed the Yang-Mills 
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FIG. 1. The Feynman rules for a pure Yang-Mills 
theory (all momenta flow into the vertices), 
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theory formally, that is without regard to the nec- 
essary renormalization procedure. There exist 
two ways of regularizing these theories. The most 
elegant method is that of dimensional regulariza- 
tion as discussed by ’t Hooft and Veltman.** Al- 
ternatively, one can add gauge-invariant higher- 
derivative terms to the Lagrangian, following 
Slavnov™* and Lee and Zinn-Justin.?5 Both of these 
methods have the necessary virtue of maintaining 
explicit gauge invariance. 

With such a regularization the only primitively 
divergent amplitudes are the vector-meson and 
ghost two- and three-point functions and the vector- 
meson four-point function. The necessary subtrac- 
tions are severely restricted by the Ward identities, 
which render all divergences logarithmic and re- 
late the various divergent amplitudes. This is 
merely a reflection of the gauge invariance which 
limits the renormalized couplings to those dis- 
played in the Lagrangian (2.5).%° 

The theory is then determined by specifying the 
subtraction constants at some convenient subtrac- 
tion point. It is, of course, impossible to subtract 
at zero four-momentum, due to the infrared sin- 
gularities. One therefore performs the subtraction 
at an (arbitrary) Euclidean point p? =- u?. 

We define the vector wave-function renormaliza- 
tion constant Z, in terms of the unrenormalized 
transverse vector propagator: 


i { kuk, 
Di (RN, =— a „+ Puka) 6.5; (2.6) 


k?=-p2 H 
the ghost wave-function renormalization constant 
Z, in terms of the unrenormalized ghost propaga- 
tor: 


ae 
Git (=u?) == 75 23593 (2.7) 


the vector charge renormalization constant Z, in 
terms of the unrenormalized vector three-point 
vertex: 


Tin’ (2,95 7)uvrl on 92 a= ZT Rape (Ps spr 


(2.8) 


and the ghost-vector charge renormalization con- 
stant Z, in terms of the unrenormalized ghost- 
vector three-point vertex: 


Ar Vee AT 


ri (ba, r)a] 
p?= 


(2.9) 


The Ward identities?t:3? ensure that these con- 
stants are related by 


A EENAA (2.10) 


and that the longitudinal part of the inverse vector 
propagator is unrenormalized: 


D i = 2S ERs (2.11) 


where a, is the unrenormalized gauge parameter. 
The renormalized Green’s functions are then de- 
fined by scaling the fields according to 


(B4), =Z, (Bia, 


2 (2.12) 
(9° ),= Zango P ’ 
and defining the renormalized charge to be 
BZ E (2.13) 


and the renormalized gauge parameter to be 
a,=Z, æ. (2.14) 


The renormalized one-particle irreducible (1PI) 
Green’s functions Rees -+P.) [T is the 
inverse propagator |, 


Dag (Pee) 


7.02 Dy ...p (Pye P,) (2.15) 


unrenormalized 9 
are then finite functions of the renormalized 
charge g, the gauge parameter a, and the renor- 
malization point u. (When g and @ appear without 
subscripts they will refer to their renormalized 
values.) 

The choice of renormalization point, u, is arbi- 
trary. Any change in u can be compensated by a 
corresponding change of the charge, the scale of 
the fields, and the gauge parameter. The renor- 
malization-group equations reflect this fact. These 
equations are most simply derived?” by noting that 
the unrenormalized 1PI Green’s functions 
p” (A; g,) (we suppress the momenta and vector 
index labels of these 1 PI Green’s functions), when 
expressed as functions of the cutoff, the bare cou- 
pling constant g,, and the gauge parameter a,, are 
independent of u: 


2 
du 


Using Eq. (2.15) and the chain differentiation rule 
we have 


U— P(A, g,, @)=0. (2.16) 


use alg, a) ae ny( g, a+ 5(g, a= | 


xP" 2, a, u)=0, (2.17) 


where 
ð 
B(g, a) = wo (2.18) 
H guy A fixed 
d1nZ 
y(g, a)= 3u = , (2.19) 
0 u Zy 1%y A fixed 


| co 


8(g,0)= p=" (2.20) 


Ə u Ey Ly ,A fixed i 
The fact that 8, y, 6 are finite functions of g and a 
is an immediate consequence of (2.17) and the fact 
that r” is finite when expressed in terms of the 
renormalized parameters. 

The last term in the renormalization-group equa- 
tion is peculiar to gauge theories. As a conse- 
quence of Eq. (2.14) we have that 


Q ƏZ 
5(g, a)=— oO eats 
Z; OU fy Ou A fixed 


=—2ay(g,a) . (2.21) 


The renormalization-group equations assume a 
particularly simple form in Landau gauge where 
a,,=a@=0, for in that gauge 6=0. This is a reflec- 
tion of the fact that the longitudinal part of the vec- 
tor propagator vanishes in this gauge, anda 
change in renormalization point does not change 
the gauge parameter. In the following we shall 
often restrict ourselves to this gauge. Ultimately 
we are interested in the physical consequences of 
these theories which are determined by gauge- 
invariant amplitudes. For such amplitudes the 
change in gauge parameter can be reabsorbed by 
a change in coupling and scale of fields. This is 
explained in some detail in Appendix A, where we 
also show that the lowest-order terms (of order 
g*) in B(g,q@) are independent of a. 

The utility of the renormalization-group equations 
is that they determine the change in the Green’s 
functions as we scale all momenta uniformly. Con- 
sider the 1PI amplitudes 


[Opi Abes- e3 PD) 


(in Landau gauge), where p; is some set of nonex- 
ceptional Euclidean momenta, and à is a nonvan- 
ishing parameter. Pure dimensional analysis im- 
plies that 


n Ş À À 
y" nDis eae »APni 8) u)= uÍ npin( MPa, Pye Abs 


I m 
(2.22) 
so that (2.17) can be rewritten as 
PÈ -eei tnie) 
an ag 
xP rAp, «-APpi Zw) =O. (2.23) 


The general solution of (2.23) is most transpar- 
ently expressed in terms of the effective coupling 
constant git, g): 


t=1na 
(2.24) 


“at, g) =6(g), g(0, g) =g. 
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This function is given implicitly by 


€ dx 


; ata) (2.25) 
and satisfies 
ð ð eee 
eDi ale, g)=0. (2.26) 


The physical meaning of the effective coupling 
constant o(t, g) is that it equals the renormalized 
coupling constant defined by performing the sub- 
tractions indicated in (2.6)—(2.9) at the Euclidean 
point 


D=- pR pe. 


It is expressed in terms of the renormalized 
coupling constant g (which was determined by sub- 
tracting at p? =— u?). The renormalization-group 
equation (2.24) then determines the effect on g of a 
change in the subtraction point. 

In terms of g we have 


POD SIASI, Pus 8) u) =T p, aoe ‘pe glt, g), u)atT” 


xexp (-n f dx y( g(x, 2) 


(2.27) 


Of particular interest is the large-) limit of this 
solution, for it determines the ultraviolet behavior 
of the Green’s functions even in the presence of 
mass terms (which here must arise from symme- 
try breaking). This limit will in turn be controlled 
by the large-¢ behavior of the effective coupling 
constant g(t, g). If Eq. (2.4) admits a solution such 
that 


lim g(t, g)=Ze, (2.28) 


t — 9% 


then we say that g. is an ultraviolet-stable fixed 
point. The asymptotic behavior of r”) is then con- 
trolled by gə according to 

Dp, ++ Pa 2 H) 


Ee i a 


xexp(-n [Pr ergo] dx) (2.29) 


a r™(p,, oa 


so that y( g.) is the anomalous dimension of the 
field. 

The fixed points of the renormalization group 
are determined by the zeros of (g); i.e., B(g.) 
=0. However, not all such zeros are UV-stable. 
Thus if £ has a simple zero at g. this will be UV- 
stable if and only if 


Bl go) =0; T Bl go) <0. (2.30) 
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A zero of § at which d8/dg>0 is said to be an in- 
frared- (IR) stable fixed point, since g approaches 
such a fixed point when A — 0 (¢— -%). 

A theory is said to be asymptotically free if g% 
vanishes. In that case y(g.)=0 and the Green’ s 
functions can be expanding for large Euclidean mo- 
menta, in a asymptotic series in g(t) [ plus a 
modification due to the integral in (2.29)]. Since 
in all theories 6(0)=0, the origin of coupling-con- 
stant space is either UV- or IR-stable. It has 
recently been proved*® that no renormalizable field 
theory without non-Abelian gauge fields can be 
asymptotically free. 

Asymptotically free field theories are clearly of 
great theoretical interest. They provide one with 
models in which the asymptotic behavior of am- 
plitudes is calculable by ordinary perturbation the- 
ory. In addition there appears to be evidence, ex- 
perimental and theoretical, that such theories are 
required to explain deep-inelastic scattering. The 
phenomenon of scaling predicted by Bjorken’ is, 
up to logarithmic corrections, a true asymptotic 
feature of asymptotically free theories. Further- 
more, it now appears that Bjorken scaling can 
only occur if the strong interactions are asymp- 
totically free.”* The fact that the only theories 
that can be asymptotically free are those involving 
non-Abelian gauge fields and that, as we shall see 
in the following, asymptotically free gauge theo- 
ries can be constructed, is a strong argument for 
a gauge theory of the strong interactions. 


Ill. CALCULATION OF THE RENORMALIZATION-GROUP 
PARAMETERS 


We shall now proceed to calculate the renormal- 
ization-group parameters ß and y to lowest non- 
trivial order in perturbation theory, that is to or- 
der g? and g”, respectively. To calculate these 
functions one must calculate the renormalization 
constants Z, and Z, to order g? and use Eqs. (2.18) 
and (2.19). These constants can only depend on the 
renormalization point u via the ratio A/p, where 
A is the ultraviolet cutoff. It is therefore sufficient 
to calculate the logarithmically divergent terms 
(of order g?) in Z, and Z,. It then follows that 


E ə Zay 
Blg, a)==8 Sn ( Zz)? (3.1) 
1 ə 
vy(g, ay) SINA InZ, ` (3.2) 


(We denote the anomalous dimension of the vector 
mesons by yy.) The calculations are thus greatly 
simplified. In particular it is not necessary to 
specify the regularization method employed. 

The wave-function renormalization constant is 
determined by the vector-meson self-energy 


| co 


(c) a~n 


FIG. 2. The vector-meson self-energy graphs. 


graphs shown in Fig. 2. We calculate from these 
graphs and (2.6) that 


2 

Z,21+-2 E- a)C,(G)InA, (3.3) 
16r 

where C,(G) is the value of the quadratic Casimir 

operator for the adjoint representation of the gauge 

group G. Namely, 


2 Caca Coca = Cz(G)5 a0 3 (3.4) 
a 


in the case of SU(N), C,(SU(N))=N. 

Similarly the charge renormalization constant is 
determined by the Feynman graphs shown in Fig. 
3 and yields 


2 
zr (G52) C,(G)lnA . (3.5) 
A useful check on the above calculations is pro- 
vided by the Ward identity, Eq. (2.10), which re- 
lates the vector-meson and the ghost-renormaliza- 
tion constants. The latter are much easier to cal- 
culate. The relevant graphs are shown in Figs. 4 


(a) 


(b) +2 crossed graphs 
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FIG. 3. The trilinear vector-meson vertex correc- 
tions, 


(c) crossed graph 
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FIG. 4. The ghost self-energy graph. 


and 5 and result in 


E E A E 

ART E 2)c,(a)ima, (3.6) 
5 g’ 

Z,=1+ 1672 (7 WEG) Ind. (3.7) 


The calculation of 6 and yy is then mere arith- 
metic. The result is 


Blg, a) =- £5 CC) +02") (3.8) 
ylas, a) = Es (8 - a)C,(6). (3.9) 


It is obvious from the above that for these gauge 
theories the origin is UV-stable. If g° is small 
enough then the solution to (2.24) will approach 
zero, as 1/t for large t. We shall explore the 
consequences of this following the incorporation of 
fermions. We cannot provide a deeper understand- 
ing of why non-Abelian gauge theories are so dif- 
ferent from all other field theories in this respect. 
Perhaps the serious infrared singularities of Yang- 
Mills theories are reflected in the IR instability, 
and thus the UV stability, of the origin. 


IV. INCORPORATION OF FERMIONS 


We now consider a gauge theory which includes 
matter fields. Fermions are easily incorporated 
without introducing any new coupling constants, 
and without destroying the asymptotic freedom. 
We add fermions by adding to the Lagrangian the 
term 


Ly = (id -M -go + B*)y , (4.1) 


where the o° are the matrices of the representation 
R of the gauge group G according to which y trans- 


(a) NZ 


FIG. 5. The ghost~ghost—vector-meson vertex cor- 
rections. 
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forms. The fermions may have nonvanishing, but 
symmetric, masses since we are not considering 
here chiral gauge groups. 

In the presence of a fermion mass the renormal- 
ization-group equations are no longer valid. In- 
stead one can derive the Callan-Symanzik equa- 
tions”! which contain, in addition to the renor- 
malization-group operator, an inhomogeneous 
term arising from mass insertions. The renor- 
malization-group parameters are unaffected by 
these mass terms. For large Euclidean momenta 
the Callan-Symanzik equations reduce to the re- 
normalization-group equations. This is a conse- 
quence, to any finite order in perturbation theory, 
of Weinberg’s theorem. More precisely consider 
the large-) limit of a Green’s function I”(p,) for 
nonexceptional Euclidean momenta. The leading 
power in A of r” to any finite order in perturba- 
tion theory is denoted by T, (àp;) and satisfies 
the renormalization-group equation [e.g., Eq. 
(2.23)]. We should emphasize that an important 
assumption is being made here, namely, that the 
leading-power behavior of perturbation theory is 
identical with that of the actual solution. 

The effect of the fermions on £ and y is easily 
calculated. The fermions contribute to Z, via the 
graph shown in Fig. 6 and to Z, via the graph 
shown in Fig. 7. These yield the contributions 


2 
Zi =Z] =- 25% T(R) In, (4.2) 
where T(R) is defined by 
Tr(o%0")=T(R)5,,- (4.3) 


We note the elementary identity between the value 
of the Casimir operator for the representation R 
and T(R) given by 


7T(R)=d(R)C,(R) , (4.4) 


where d(R) is the dimension of the representation 
R and v the dimension (number of generators) of 

the group. For example, in the case of SU(N) we 
have for the vector representation (N) 


T((N)) =3, 


Nea 
2N ’ 


C,((N)) = 


whereas for the adjoint representation T(ADJ)=N 


=C,(G). 


FIG. 6. The contribution of the fermions to the vector- 
meson self-energy. 
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FIG. 7. The contribution of the fermions to the trilin- 
ear vector-meson coupling. 


The equality of Z¥ and Z? is an immediate con- 
sequence of the Ward identity, Eq. (2.10), and the 
fact that Z, and Z, receive no contributions from 
fermion loops to lowest order. The resulting con- 
tributions to (g) and y( g) are 


Bla) =+ STR), (4.5) 


T E , 8 ap) (4.6) 
4 167? 3 " 


These contributions are opposite in sign to those 
arising from the vector-meson loops. Thus the 
fermions tend to destabilize the origin. However, 
there is room to spare. As long as 


T(R) < 4# C,(G) (4.7) 


the theory will be asymptotically free. This re- 
quirement is not very restrictive as to the number 
of fermions allowed in the theory. For example if 
= the gauge group is SU(3) (C,=3), one can accom- 
modate as many as 16 triplets or 2 octets of fer- 
mions without losing asymptotic freedom. 

It is therefore possible to construct a large class 
of theories in which the renormalization-group 
equations for the effective coupling constant g(t, g) 
take the form 


d — ss dle es = 
Ape =2°8(g)=-bog* tbig > tee, (4.8) 
where 
1 
bo = gals Ca(G) -3 T(R)]. (4.9) 


The structure of this equation assures us that as 
long as g(0, g) is small enough then for large ¢ we 
will have 


gl, g) = bat T+ biba t 7 Int+ O(1/t?) . 
t- 


(4.10) 


This will be the true asymptotic behavior of the 
effective coupling constant as long as g(0, g)=g is 
in the domain of attraction of the fixed point g = 0.°° 
The size of this domain is determined by the value 
of the first zero of (g), namely, if 


B(g)=0 for g=g, , (4.11) 


then the domain of attraction of the origin is given 
by?’ 

O0<g*<g,’. (4.12) 
Perturbation theory of course tells us very little 
about the nonvanishing zeros of 8. If we calculate 
6 to any finite order, we might expect zeros to oc- 
cur at values 


2” OC) =O(1) 


re (4.13) 


since this is the effective expansion parameter. 
On the other hand such approximations to 8 are 
totally unreliable for these values of g. It is per- 
fectly possible that in gauge theories 6 is negative 
semidefinite and all values of the coupling constant 
are in the domain of attraction of the origin. 
Clearly for such theories to describe the strong 
interactions it is necessary that this domain be 
relatively large. It would therefore be useful to 
know the value of 6,, which we have not calculated. 

The value of b, would be interesting in another 
context. Since it is possible, by including the req- 
uisite number of fermions, to render b) very 
small, one might hope to construct models for 
which g, is very small. For example if the gauge 
group is SU(3) and we have 16 or 17 triplets of 
fermions, then b, will equal 1/247? or ~1/247’, 
respectively. This value is suppressed by a factor 
of roughly 30 compared to the “natural” scale of 
bọ- Therefore unless there are similar cancella- 
tions in the calculation of b,, we would expect g,” 
to be rather small, and calculable to a good ap- 
proximation from the two-loop expression for £. 

If this is so one could construct models which have 
UV-stable fixed points at zero or g, and IR-stable 
fixed points at g, or zero, respectively. These 
would provide interesting theoretical models in 
which both the ultraviolet and the infrared asymp- 
totic behavior would be calculable. 

The physical consequences of asymptotic free- 
dom will be explored in the following section and 
in a subsequent paper. It is clear, from the discus- 
sion in Sec. II, that the ultraviolet asymptotic be- 
havior of all Green’s functions can be calculated. 
Thus the z-vector-meson 1PI Green’s functions 
will behave, for large Euclidean momenta, accord- 
ing to Eq. (2.29). If we define (we shall work in the 
Landau gauge) 


yyle)=cog +g? ters, (4.14) 
where 
| r 

CaS TE C(G)-ŻT(R)], (4.15) 


then 


| œ 


rrp, oe Abang Le) Ai A47" (Ina) 770 do I, 


i 1 Ininà 
x [ri (Dosser pe O DAT U peba b max? to ie 
(0) 


where 


l g g 
I,=exp(-n f dl ry x 8)) = CoB) m~ const+O 
0 


A> 0 


In particular, the transverse part of the vector- 
meson propagator behaves like 


Ruky P 
Drogo anon ti ann, 


(4.18) 


To calculate the asymptotic behavior of Green’s 
functions involving fermions we require the anom- 
alous dimension of the fermion field yp. This is 
readily calculated (except for group-theoretical 
factors it is the same as in QED): 


yr(g a) =g fott" 


g’ 
“Fee LCA) tor? : (4.19) 
Thus the fermion propagator behaves for large mo- 
menta according to 


Sp(p) ~ gin pry! bo (4.20) 


b —> = 00 

It is an important feature of asymptotically free 
theories that the Green’s functions are not strictly 
given by their free-field-theory expression for 
infinite momenta, due to the presence of the loga- 
rithmic term in (4.16). These arise because the 
anomalous dimensions vanish like g°, and g° van- 
ishes only logarithmically as A— œ. This, of 
course, is a consequence of the fact that B(g*) has 
a double zero at g?=0. These logarithmic correc- 
tions to free-field behavior are calculable and in- 
dependent of the coupling constants. They can be 
computed by calculating the relevant anomalous 
dimension to second order. The above Green’s 
functions are of little physical interest due to their 
gauge dependence. However, as we shall see 
shortly, the physically interesting Green’s func- 
tions of gauge-independent operators (currents in 
particular) exhibit essentially the same asymptotic 
properties. 

We emphasize that the derivation of the Green’s 
functions in asymptotically free theories does not 
require the existence of a convergent perturbation 
theory. It is sufficient to assume that perturba- 
tion theory yields an asymptotic expansion, for 
small coupling constant, of the relevant ampli- 
tudes. In that case the renormalization-group 
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ae (4.16) 


: ) l (4.17) 


a 


equations provide, for asymptotically free theo- 
ries, a true asymptotic expansion of Green’ s func- 
tions for large Euclidean momenta. 


V. MODELS AND APPLICATIONS 


In order to construct realistic non-Abelian gauge 
models of the strong interactions one must con- 
front the issue of symmetry breaking. The stan- 
dard means of breaking the gauge symmetry is to 
introduce scalar mesons explicitly into the La- 
grangian. As will be explained in Sec. VI it is very 
difficult to preserve asymptotic freedom while in- 
corporating scalar mesons, and perhaps impos- 
sible to include enough scalars to completely break 
the gauge symmetry. This is not too disappointing 
since the explicit introduction of scalar mesons, 
whose only role is to break the symmetry, is not 
very pleasing. An alternative is that the gauge 
symmetry is dynamically broken. In other words, 
a composite Goldstone boson is formed and is 
eliminated by the usual Higgs mechanism. *° 
Another possibility is that the gauge symmetry ts 
exact. At first sight this would appear ridiculous 
since it would imply the existence of massless, 
strongly coupled vector mesons. However, in 
asymptotically free theories these naive expecta- 
tions might be wrong. There may be little connec- 
tion between the “free” Lagrangian and the spec- 
trum of states. 

The possibility of dynamically induced spontane- 
ous symmetry breaking has been considered by 
many authors,*”” although no realistic model has 
been constructed which does not involve fundamen- 
tal scalar fields. Of particular interest is the 
work of Coleman and Weinberg.* They show that 
theories involving massless particles often become 
unstable, due to infrared singularities, and ex- 
hibit spontaneous symmetry breaking. The infra- 
red singularities of a Yang-Mills theory are par- 
ticularly severe. Furthermore, zero coupling, 
for such theories, is an ultraviolet-stable fixed 
point and therefore infrared-unstable. This means 
that (neglecting masses) as the momenta decrease 
the effective coupling constant increases. Pertur- 
bation theory is therefore totally unreliable inso- 
far as the small-momentum behavior of an asymp- 
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totically free theory is concerned. The same re- 
normalization group techniques allow one? to dis- 
cuss the small- or large-classical-field behavior 
of the “potential.” The infrared instability of the 
origin indicates the unreliability of the classical 
(free) approximation to this potential. Thus wheth- 
er or not the theory exhibits symmetry breaking is 
a difficult dynamical question, requiring nonper- 
turbative calculations. 

If the gauge symmetry is broken by a dynamically 
induced Higgs mechanism then the vector mesons 
will acquire masses (say 1—3 BeV) and the color 
degeneracy will be split. In that case one would 
still be faced with the fact that there is no experi- 
mental evidence for the existence of such neutral 
vector mesons, colored hadrons and especially 
quarks. The proponents of “red, white, and blue” 
quarks” as a mathematical abstraction argue that 
the color SU(3) group should be exact, and that all 
noncolor singlets should be suppressed completely. 
One clearly requires a dynamical explanation of 
such a miracle. It might very well be the violent 
infrared singularities of an asymptotically free 
gauge theory provide the requisite dynamical 
mechanism. 

To illustrate some of the strange things that 
could occur in the region of small momenta con- 
sider a pure Yang-Mills theory (no fermions). 
Then the renormalization group equations, Eq. 
(2.23), as well as the general solution, Eq. (2.27), 
are exact. One can use them to discuss the on- 
mass shell, or infrared behavior of the theory by 
letting à— 0, ż=lnìà— ~œ. The effective coupling 
constant g, which controls the dynamics in this 
region is given [in terms of g(t=0, ¢)=g] as usual 
by 


g(t) dx 
J OME (5.1) 
The behavior of g as ¢--—% will depend on the ac- 
tual form of (g). We can distinguish two cases: 

A. B(g) vanishes at g=g,;<~. In this case g(t) 
will approach g, as t+--©, its rate of approach 
will depend on the nature of the zero. This is the 
simplest case to envisage, the Green’s functions 
will scale according to Eq. (2.29) with some anom- 
alous dimension. 

B. B@) is always negative. Here we must fur- 
ther specify whether the integral 


” dx 
T J T (5.2) 
is finite or not. 

If T=œ, then g(t) approaches infinite values as 
t--—©, Indeed if B(x)=x°, œ< 1, for large x then 
glt) ~ (=t) C-% for large t. 

If on the other hand T is finite then the effective 


coupling constant diverges for finite momenta. In 
fact, g(T,g)=%. 

The infrared behavior of Green’s functions in 
this case is determined then by the strong-coupling 
limit of the theory. It may very well be that this 
infrared behavior is such so as to suppress all but 
color singlet states, and that the colored gauge 
mesons as well as the quarks could be “seen” in 
the large-Euclidean-momentum region but never 
be produced as real asymptotic states. This is an 
exciting possibility which requires further exami- 
nation. ** 

In any case it might be valid to assume that what- 
ever happens to the theory for small momenta 
does not affect the ultraviolet behavior. We shall 
therefore construct models and calculate quantities 
of physical interest leaving the problem of symme- 
try breaking (or the lack of it) to further work. 

Since we are not to worry about symmetry break- 
ing, our models need not include scalar mesons. 
We therefore have only to specify the strong-in- 
teraction gauge group, G, and the fermion repre- 
sentation R. We would like, of course, to pre- 
serve Gell-Mann’s (approximate) SU(3)xSU(3) sym- 
metry. This is simply achieved by taking the 
gauge group to commute with the ordinary SU(3) 
xSU(3) generators, and having the fermions belong 
to a representation of SU(3)xSU(3)xG.* We shall 
take the fermions to be ordinary triplet quarks. 
The fermions can then be represented by a matrix 
of spinor fields 


W=f 0 7 7 4, (5.3) 


Cr Ny Àr 


The generators of SU(3)}xSU(3) transform the col- 
umns of this matrix, whereas the generators of G 
transform its rows. The fermions thus transform 
as an L-dimensional representation (R) of the 
gauge group G. 

In such a scheme the vector mesons associated 
with the generators of G are neutral with respect 
to SU(3)xSU(3). The labels 1,2,..., L which dis- 
tinguish the different quark triplets can be thought 
of as colors, so that the strong interactions are 
mediated by colored gauge mesons. Colored 
quarks have been considered before* for other 
reasons, and there is some evidence that three 
colors would be welcome. We shall therefore con- 
sider a model in which the strong gauge group is 
SU(3) and the fermions are color triplets (ZL =3) 
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(although as far as asymptotic freedom is con- 
cerned any group will do). 

In this model © transforms under an infinitesimal 
gauge transformation according to 


W(x) ~ U(x) +i la U(x) (5.4) 


and under an ordinary SU(3)xSU(3) transformation 
according to 


W(x) —~ U(x) +i Elx) E, (5.5) 


[the à? are the usual SU(3) matrices]. Our Lagran- 
gian is 


L=Ly+Tr{ pis - gB, X) - PYM}, (5.6) 


where Ly is given by Eq. (2.5) and M is the fermion 
mass matrix. This model is asymptotically free. 
The numerical value of 6(g) is 


Bla) = -2g 


Aes (5.7) 


The ordinary SU(3)xSU(3) vector, V, and axial- 
vector, Aj, currents are given by 


Ve =Tr{ y ya}, 


AS, =Tr{ Ey Ysp" } - 


Let us now discuss deep-inelastic scattering in 
such an asymptotically free theory. Due to the fact 
that the vector mesons are neutral and that the the- 
ory is asymptotically free one can derive all the 
sum rules and relations previously derived in the 
parton or light-cone models. In addition there will 
be logarithmic deviations from Bjorken scaling 
which can be calculated. The full analysis will be 
given ina subsequent paper *°; here we shall give a 
simplified discussion. 

The structure functions of deep-inelastic scatter- 
ing are Fourier transforms of the product of elec- 
tromagnetic currents. In the scaling region one 
probes this product for lightlike separation of the 
currents. To discuss this one employs Wilson’s 
operator-product expansion: 


Yo 
2 
g? 

Da 


(5.8) 


J (xJ (=x) ~ 5 Cx? g) yte «extn O s: „„ (0) 


x ro n 


(5.9) 


[where we have suppressed the vector and SU(3) 
labels on the currents as well as the tensor and 


SU(3) structure of the operator-product expansion]. 


The dominant operators in the scaling region are 
those of twist (equal to dimension minus spin) 2. 
These are denoted by O™, The c-number function, 
CXx?, g), contains the light-cone singularity and 
controls the asymptotic behavior of the structure 


functions. In fact, the moments of the scaling 
structure functions measure the Fourier transform 
of c: 


1 
f dxx"™F laa) o Pla, gX2|O™|n) 
0 


-q -> 00 
(5.10) 


(where xis the standard Bjorken variable, x 
=~q?/2v, and F; stands for vW, or xW,). Of 
course the moment is also proportional to the nu- 
cleon matrix element of the operator O}; however, 
the q? dependence is contained in C'"\q?, g). 

One can apply renormalization- group techniques 
to Wilson’s expansion to derive an equation relat- 
ing the dependence of C'” on x? and g.?!~*3 If the 
operator O' is multiplicatively renormalizable 
(renormalization constant Z,), then one derives 


ə ə zoa 
atel 2y 7(g) - yala) Cc (Ss) =0, 


(5.11) 


where y, is the anomalous dimension of the “cur- 
rent” J(x), and y, the anomalous dimension of 
ol”. 


Q 
Ya (D) =u nZ, (5.12) 


OM Gy Ou A fixed 
In the case under consideration the currents 
(linear combinations of V} and Aj) are conserved 
or partially conserved. They therefore have van- 
ishing anomalous dimension, so that the solution 
of Eq. (5.11) is 


2 =< 
(4, «) = P, a(t, 2)) 


\ 


x exp(- f. (Bx, eae), (5.13) 


In our asymptotically free models, g— 0 and y,( g) 
vanishes as t- ©. However, y,(¢) does not vanish 
fast enough to render the integral in Eq. (5.13) 
convergent. Let 


Vn 2°) =Vne? + Ol g*); (5.14) 


~1 


then since g?°— btt, we have that 


; 2 v 
o(a) ~ CP (1, 0)0ng?) %2% 
xexp(- f [ra (Ee - v,g?] ds). 

O 


(5.15) 


Therefore the logarithmic deviations from Bjorken 
scaling are obtainable from the lowest-order cal- 
culation of y, (¢7). 

In general a given operator o™ will not be mul- 
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tiplicatively renormalizable. In particular if there 
exists more than one operator with the same spin, 
quantum numbers, and physical dimension then 
one must take linear combinations of these to ob- 
tain operators with definite dimensions. This mix- 
ing occurs for the twist-2, SU(3)xSU(3) singlet 
operators in our theory, since there exist more 
than one such operator. Thus, for example, for 
spin 2 both 


OF = Trivyy (3, -2gB; ha)bt 


+(u—v)-— (trace terms) (5.16) 


( 
OT Gt) 


Ey Triat F(x, 


where V, =, +ig BGA" is the covariant derivative. 
We define Cota to be an 1PI Green’s function with 
the insertion of the operator One. u(x). This am- 
plitude will satisfy a renormalization-group equa- 
tion. In particular if O'”* is inserted into the fer- 
mion two-point function we obtain 


ugg elegy rila- 2y (2) rG ata =0, 


(5.18) 


where yp is the anomalous dimension of the fer- 
mion field. Since £ is of order g? we may calcu- 
late the combination y$ — 2yp to order g” by evalu- 
ating the logarithmically divergent contributions 
to Pa to order g°. These are given by the Feyn- 
man “graphs of Fig. 8 and yield 


yele) =yig? t... 


2 2 n 
= 23 C,(R) É > 


1 
nin+ 1)" aa 


(5.19) 


lin the three-triplet model C,(R)=4+]. As expected 
the anomalous dimensions are independent of a and 
y. Having performed this calculation we can now 
compute the scaling behavior of the nonsinglet 
pieces of the deep-inelastic structure functions ac- 
cording to Eq. (5.15). If we denote by N(x, q7) one 
of these structure functions, say Fproton — ppeutron 


FIG, 8, The lowest-order correction to re) ae 


Mites vie” peas 


| co 


and 


(2) 1 
O Sik. = A age, 


contribute to the operar product expansion a 
=ð, B} -Ə Bh -ig fave BIBS). This mixing and the 
Saiculation of the resulting anomalous dimensions 
will be discussed in a forthcoming publication.*° 

The nonsinglet SU(3)xSU(3) twist-2 operators 
however are uniquely given in terms of the fermion 
fields, and have definite dimensions. Let us denote 
these operators by O14. (x,7), where r (=+1) 
denotes the chirality of tie. operator. They are 
given by 


"Vy (1+ry,5)%(x)} - (trace terms), (5.17) 
E 
then 
fiar N(x, q) œ~  const(Ing?)7%n iof :)|, 
g2~> wo Ing? 
(5.20) 
where 
3C,(R) 2 ‘ | 
C ee a eco ee oe = 
A, fe o arin A kj. 
(5.21) 


We note that the power of the log is independent of 
the physical coupling constant and is determined 
solely by the gauge group (G) and the fermion rep- 
resentation (R). For the three-triplet SU(3) model 
A =}, 4,=23, and for large n, A, =21n(z). 

In a forthcoming publication we shall analyze 
lepton-hadron scattering in great detail; meanwhile 
a few comments are in order. 

a. The approach to asymptopia in these theories 
is logarithmic. All asymptotic relations will be 
corrected up to terms of order 


i 2 \=1 
b,~* (nts) ; 
The rate of approach to asymptopia depends on the 
unknown low-momentum behavior of the theory. 

b. Bjorken scaling is violated by logarithmic 
terms, as in Eq. (5.20). These logarithmic correc- 
tions are given by second-order perturbation-the- 
ory calculations, and depend only on the gauge 
group and the fermion representation. The only 
case where such logarithms are absent is when 
the relevant operator in the Wilson expansion has 
zero anomalous dimension. This is the case for 
vector and axial-vector currents as well as the 
energy-momentum tensor. Thus the first moment 
of the isoscalar component of the structure func- 
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tions does scale: 


1 E 1 ) 
I=0 Syd 
J dx F1™ (x; q ae const +0( i i 
(5.22) 


whereas the higher moments decrease logarithmi- 
cally as in Eq. (5.20). Similarly the total neutrino- 
hadron cross section will scale, since it is 

given by the same (spin-2) moment. 

It is an open question whether such a picture is 
consistent with experiment. The deviations of the 
asymptotic forms of the moments from exact 
scaling are quite small, and one would need rather 
large variations of q? to see any marked change. 
Verification of logarithmic deviations, as well as 
the numerical values of a,, would be strong evi- 
dence for the above class of theories. 

c. Sum rules and other relations between the 
various structure functions measured in electron 
or neutrino scattering, which previously were de- 
rived from current algebra, the parton model, or 
a naive light-cone expansion, are all true theorems 
in our models. This is because, as is seen in Eq. 
(5.15), the SU(3)xSU(3) and tensor structure of the 
Wilson expansion will be (for large q°) that of the 
free-quark model—the coupling constant effective- 
ly vanishing in this region. There is no point in 
listing these predictions here—they have been re- 
viewed by many authors.**:* However, we note 
that these relations are approached logarithmical- 
ly. Thus, for example, since the charged consti- 
tuents have spin $ (Ref. 16) 


1 9, n 
SF (x, q°)x dx 1 
[OF p(x, q?)x"dx q? > =o AA Gane Ms 


(5.23) 
where 


1 
~~ —. 
Goes Ing” 


f (lng?) 


d. Electron-positron annihilation at large ener- 
gies is controlled by the identity operator in the 
Wilson expansion of two electromagnetic currents. 
Since these all have cannonical dimensions the 
cross section scales, and the coefficient is deter- 
mined by the free-field limit (s = center-of-mass 
energy squared) 


O 


ete- —> hadrons _. 2 S J 
= . ^| 1 +O ln— . 
Opie ea 2,9 ( z) 


Finally let us discuss the incorporation of weak 
and electromagnetic interactions into our models. 
After all the recent revival of gauge theories was 
for the purpose of constructing a unified and re- 
normalizable theory of the weak interactions.*® Is 
there any problem in combining these theories? 


(5.24) 


Since in our models the strong gauge group com- 
mutes with SU(3)xSU(3) one can easily incorporate 
the weak plus electromagnetic interactions accord- 
ing to any one of the various schemes proposed 
recently. In fact, as far as the weak interactions 
are concerned there is no difference between our 
models and an Abelian vector-gluon interaction. 

One would have to go to extraordinarily high 
energies to ascertain experimentally whether the 
weak or electromagnetic interactions are asymp- 
totically free [e?In(q?/M,”)=1]. If however we as- 
sume that these theories are asymptotically free 
(in which case one is restricted to semisimple 
gauge groups and one must worry about symmetry 
breaking as before) then the Baker-Johnson-Adler 
approach to QED’° would be unnecessary. The 
ultraviolet behavior would be controlled by the 
fixed point at zero coupling. 


VI. THE INCORPORATION OF SCALAR PARTICLES 


We shall now consider non-Abelian gauge theo- 
ries which include scalar fields. There are a 
variety of reasons for considering such theories. 
One might want to incorporate spin-zero funda- 
mental fields as such, or one might want to employ 
such fields to provide masses for the vector me- 
sons by means of the Higgs mechanism. Indeed 
the only known way of breaking the gauge symmetry 
spontaneously is via scalar mesons which develop 
nonvanishing vacuum expectation values. We have 
explored whether one can incorporate a sufficient 
number of scalar mesons to break the gauge sym- 
metry completely, or at least retain only an 
Abelian gauge group. In both these cases one 
would thereby have a theory which is both asymp- 
totically free and which possesses an S matrix in 
perturbation theory (and not just off-shell Green’s 
functions). 

It is well known from the proofs of the renormal- 
ization of Higgs theories* that their ultraviolet 
behavior is tha. of the underlying symmetric theo- 
ry. In other words the symmetry breaking is a 
“generalized mass term” and does not affect the 
asymptotic (Euclidean) behavior of the theory. 

The argument that this is so is implicit in the proof 
of renormalizability of Higgs theories,’ where 

one shows that the subtractions which render the 
symmetric theory finite are sufficient to make the 
asymmetric theory finite. 

Let r” p, ..+,P,3v) be any Green’s function of 
the asymmetric theory, where v is the vacuum ex- 
pectation value of a scalar field which breaks the 
symmetry. I” may be represented as a function- 
al derivative (n times), with respect to the fields, 
of the generating functional of 1PI’s. It is an im- 
mediate consequence of this definition that 
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rp, sete Dr) 
o j 
DA PrN p,, tc hgD Os rearg Oe 0), 
o J 
(6.1) 


where on the right-hand side of (6.1) the r‘"*”) are 
the Green’s functions of the symmetric theory with 
j insertions of the scalar field carrying zero mo- 
mentum. Now, by simple power counting the 1PI’s 
r*+ are less and less divergent in the ultraviolet 
region as j increases, so that the ultraviolet be- 
havior of (6.1) is controlled by the first term, i.e., 
by the ultraviolet behavior of the symmetric theo- 
ry. Thus the desirable properties of the symmet- 
ric theories discussed below will remain intact in 
the presence of symmetry breaking. 

Before we plunge into the analysis an overview 
is in order. In the case of pure gauge theories as- 
ymptotic freedom provides no constraint. When 
the theories include fermions a weak constraint 
(not too many fermion multiplets) must be satis- 
fied. On the other hand the requirement of asymp- 
totic freedom will severely constrain gauge theo- 
ries involving scalar particles. This is because 
such particles will necessarily have their own self- 
couplings. One must therefore investigate the 
asymptotic freedom of these new (dimensionless) 
coupling constants. It is well known that a scalar 
field by itself, with a \¢* coupling, is not asymp- 
totically free. Therefore the only hope is that the 
gauge mesons will help render the fixed point A= 0 
ultraviolet stable. This turns out to be very dif- 
ficult to achieve—in fact, it is remarkable that it 
is possible at all. Unfortunately in order to have 
asymptotic freedom one is forced to large gauge 
groups and representations. Furthermore we have 
not been able to find any examples of models in 
which the gauge symmetry is completely broken. 
For these reasons the results presented below 
reinforce our expectation that if asymptotically 
free theories of the strong interactions are to be 
sensible, then the symmetry breaking must be dy- 
namical. 

Let us now analyze the renormalization- group 
equations for gauge theories involving scalar me- 
sons. First we note that the scalars have essen- 
tially the same effect on g as the fermions. In the 
vicinity of the origin in coupling-constant space, 
the renormalization-group equation for the gauge 
coupling g does not involve the scalar self-cou- 
plings. These would contribute terms to Eq. (4.8) 
of order g*\*, where à is some scalar self-cou- 
pling. However, unless à <O(g*) the scalar 
self-couplings are not asymptotically free (see be- 
low) and thus these corrections are negligible. 

The net effect of the scalars is to contribute to 


B(g) a term equal to z (a)of the fermion contribu- 
tion if the scalars transform according to a real 
(complex) representation of the gauge group. Thus 


4 C,(R)d(R) 
3 y 


1 (Rath) 
6 Y 


B(g) m E BS C,(G) = 


fermions 


scalars 
=bg? ters (6.2) 


[where d(R) for the scalars equals the real dimen- 
sion of the scalar representation R], There is, 
therefore, no problem in retaining asymptotic 
freedom for g, as long as the number of scalars 
is not too big. 

In addition we are now required to study the re- 
normalization-group equations for the scalar cou- 
plings. There are, in general, many such cou- 
plings. We are, of course, only interested in the 
quartic couplings. One such coupling exists for 
any representation of the gauge group, namely the 
square of the inner product (the quadratic Casimir 
operator). If the representation in question is real 
the vertex is given in Fig. 9, and if it is complex 
the vertex is given in Fig. 10. 

A general representation of the gauge group will, 
in general, possess many additional quartic invari- 
ants. For each such invariant there will be per- 
force a dimensionless coupling constant (some ex- 
amples will be given below). Therefore, in gener- 
al, the renormalization-group equations will be a 
set of coupled nonlinear differential equations. It is 
useful to consider first a simple example in which 
there is only one quartic self-coupling. This is 
the vector representation (real dimension =2N) of 
SU(N). If @ denotes the (complex) scalar field 
transforming according to this representation, the 
scalar interaction is given by 


L= |(0, — igB,) |? - 3A(H* GY. (6.3) 


The renormalization-group equation for à will 
now be 


dn bee a 
ap Pas 8): A(t=0,rA)=a, (6.4) 
where 
Or 
Peu . (6.5) 
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FIG. 9. The canonical quartic scalar coupling for real 
representations, 
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FIG, 10. The canonical quartic scalar coupling for 
complex representations. 


Once again it is sufficient to calculate the lcga- 
rithmically divergent corrections, in the one-loop 
approximation, to the quartic scalar vertex. 
These are illustrated in Fig. 10. [The graphs in 
Fig. 11(c) should be divided by 2 since it is the 
square root of the wave-function renormalization 
constant that enters the renormalization-group 
equations.] The terms illustrated in Figs. 11(a)- 
11(d) contribute to 6, terms of order X°, rg’, Ag’, 
and g*, respectively. We therefore obtain an 
equation for the effective coupling constant i: 

dr 


— =A? + Br 97+Cg*. 


it (6.6) 


The values of A, B’, and C are readily calcu- 
lated (they are gauge~independent, but it is sim- 
plest to evaluate them in Landau gauge): 


1 
A=_a(N+4), (6.7) 
,. 1 3(N?-1) 
= 1 3(N—1)(N?+2N- 2) 


Equation (6.6) is most easily analyzed by intro- 
ducing the parameter 


a =(g)7?A; 
We note that à must be of order g’, for if it were 


larger then, since A >0, the theory would not be 


asymptotically free. The coupling a satisfies 
a E E E EEEE 


a(0)=g 7°. (6.10) 


seti 
1 2A 2A 

we will be driven to the value œ =g, as t— œ. In 

other words a, is a UV-stable fixed point. Equa- 

tion (6.11) can easily be solved by quadratures: 


alt) dx t g? 
S e a E 6.15 
J Ax? + Bxt+C as (6.15) 
For large ¢ we have 
alt)— a, m tA atza b0 (6.16) 


{-% 


and since g?*=1/t, A(t) will approach zero like 1/t. 
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FIG. 11. The graphs that contribute to the 8 functions 
for the quartic scalar couplings (the directed lines re- 
fer to the complex scalar mesons). 


GO =F? (Aa? + Ba +C), (6.11) 
where 
B=B'+b,. (6.12) 


The fixed points of (6.11) will be given by the 
zeros of the quadratic form. These will be real 
and positive if 


-(4A4C)"? <B <0. (6.13) 


When this condition is satisfied (A and C are al- 
ways positive) then there will exist two fixed points 
a, and a,. If the coupling constant à = a(0)g" is 
chosen so that 


4AC)*/?] , (6.14) 


ŘS 


It is easy to see that the condition (6.13) will be 
satisfied, for SU(N), as long as N>3 [not SU(2)] 
and 6, is small enough. We can always make b, 
small by including fermion representations, whose 
only effect is to decrease b,. (We can exclude 
Yukawa couplings by the discrete symmetry @ 
+-.) In the case of SU(3) we must incorporate 
16 fermion triplets in order to make b, equal to 
—1/487 2, 

We see already in this simplest of all cases that 
the UV stability of à is extremely delicate. Even 
in the most favorable case, N>>1 and 5,*0, we 
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can accommodate at most two vector multiplets of 
scalar fields. Moreover if we do not include 
enough fermions to render b, sufficiently small we 
cannot include any scalars. Roughly speaking un- 
less b is small the gauge mesons, which are sta- 
bilizing the inherently unstable ,¢* coupling, van- 
ish too rapidly as t~ œ. 

Clearly we cannot include enough vector repre- 
sentations to break the gauge symmetry complete- 
ly. We have also investigated many other models 
with scalar particles (i.e., other groups and rep- 
resentations). We shall just give a brief summary 
of the more interesting results. 

The adjoint and the symmetric tensor represen- 


tation of SU(N), as well as the (N, Ñ) representation 


of SU(N)xSU(N), are examples of theories with two 
scalar quartic couplings. They lead to renormal- 
ization-group equations of the form (6.6), which 
can again be simplified by dividing by g? as in 
(6.11). The result has the form 


1 
a AH AZu a+ AS gap +A BE + Bat Ca, 

ra (6.17) 
za gi Aen +AagaB+Ag sh + BB+ Cg : 


The fixed points of these coupled equations are 
easy to analyze. For each value of œ there will 

be two roots of dß/dt=0. If these roots are real 
then the smaller, 8_(a), is attractive whereas the 
larger, (a), is repulsive. Similarly we define 
a_(8) and a, (8). The coupling-constant plane is 
the pictured in Fig. 12, in the case where all roots 
are real. We are clearly driven to the fixed point 
Qs, Bs satisfying 


a _(B,) = Qs 5 
B_(a,) =P; 
as long as the initial values are small enough: 


a(0) <a ,(B(0)) , 


(6.18) 


(6.19) 
8(0) <8, (a(0)). 


In practice the simplest way to search for a 
fixed point is to solve Eq. (6.18) by iteration. We 
have mapped out some coupling-constant trajec- 
tories on a computer, a sample result is given in 
Fig. 13. 

We shall now briefly review the results for the 
above-mentioned theories. The same Feynman 
graphs as in the case of the vector representation, 
Fig. 12, must be evaluated, except that there are 
now. two types of scalar couplings. In each case 
a=)/g? and B=n/g?. 

(a) The adjoint representation of SU(N). The 
scalars are Hermitian traceless NXN matrices 


| œ 


FIG. 12. An illustration of the fixed-point trajecto- 
ries in the a-8 coupling-constant plane. The cross in- 
dicates the presence of an ultraviolet-stable fixed point. 


M of fields. The quartic self-couplings are given 
by 


L=- Tr{M?} - 4 Tr{M*} (6.20) 


(ob) The symmetric tensor representation of 
SU(N). The scalars are represented by an NXN 
symmetric matrix M of (complex) fields. Under 
the group transformation they transform as M 
>~ UMU? . The quartic scalar interactions are 


L=-\Tr{MM" ly -nTr{(um'ty}. (6.21) 


(c) The (N,N) representation of SU(N). The 
scalars are represented by an NXN matrix M of 
complex fields. The quartic self-couplings are 


L=-3x(Tr{MM"}) — $nTr{(umtyt, (6.22) 


The results are shown in Table I. The main re- 
sult is that all these theories are asymptotically 
free for N large enough. The UV stability is, as 
before, delicate. One must make b, small, and 
additional symmetrically coupled scalars destroy 
the asymptotic freedom. 

Our experience in these and other cases sug- 
gests that it is the size of a representation, rather 
than its detailed nature, that determines whether 
the model will be asymptotically free or not. In- 
deed in all the cases enumerated in Table I the re- 
normalization-group equations are identical as N 
oo, If one considers larger representations the 
numberof independent coupling constants rapidly 
increases and the renormalization-group equations 
become quite complicated. For example, we have 
analyzed a theory involving an SU(N)xSU(N) gauge 
group, with the scalars transforming according to 
the (N, N)®(N, 1)@(1, Ñ) representation. In this 
model there exist five independent scalar coupling 
constants. It is asymptotically free when N= 5. 

Our original motivation for studying theories in- 
volving scalar mesons was to utilize the Higgs 
mechanism to break the gauge symmetry and pro- 
vide masses for the vector mesons. Some of the 
models considered above would appear to have the 
potential of accomplishing this aim. Take, for ex- 
ample, the gauge group to be SU(N)xSU(N), and 
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TABLE I, The renormalization-group parameters for two-coupling-constant models. The coefficients refer to Eq. 
(6.17). The BY coefficient is evaluated for B,=0; otherwise BY must be modified according to Eq. (6.12). 


Stability 
Representation Y Aba AY g Akg BY CY criterion 
SU(N): adjoint Q n aN" = 6 SN" +18 —6N 9 N=6 
Bie 2 N N? = 
2N? —18 3 
B 0 6 N? —6N oN 
2 = 2 
SU(N): a MEA 4(N +1) 6 — e “a N2=9 
symmetric 9 9 
6(N° +N —2) 3N*+12N —48 
tensor B 0 6 2N +5 are aaa 3N 
2 Deas a 
SU(N) x SU(N): oe = i : 4N 6 -I= aN ~ L N=5 
(N,N) 
_ 6(N* —1) 3N? —24 
B 0 6 2N N oN 


the scalars to transform according to the (N, N) 
representation (which is asymptotically free for 
N>5). The scalars are represented by an NXN 
matrix M of complex fields, and they transform 
according to M — GMH". If the scalar mesons 
develop a nonvanishing vacuum expectation value, 
(0|M|0) =M., the symmetry will be broken. The 
criterion that the gauge symmetry be completely 
broken and that all vector mesons acquire a mass 
is that there be no subgroup of transformations 
that leaves M, invariant.*” If M, were an arbitrary 
NXN complex matrix this criterion would be sat- 
isfied. However, M, is not arbitrary, rather it is 
determined by the form of the Lagrangian. Thus 
-in lowest order M, is determined by minimizing 
the “potential” 


—~L,=-—p?Tr{MM"}+,. Tr{(um"y} 


+n Tri(MMŻY}. (6.23) 
Now MM" is a positive Hermitian matrix and the 
potential is a function of the squares of its eigen- 
values. When one minimizes (6.23) one finds that 
these eigenvalues are all equal in magnitude so 
that M,M; is a multiple of the identity. In that 
case M,is, upto a constant, a unitary matrix. 
There exists an SU(N) subgroup of gauge transfor- 
mations, G=G, H=M,GM,™", that leaves M, in- 
variant. Such a model will therefore contain N? 
~ 1 massive and N? - 1 massless vector mesons. 
In a similar fashion the other two asymptotically 
free models described above contain massless 
vector mesons even when the scalars have a non- 
vanishing vacuum expectation value. Basically the 
problem is that asymptotic freedom requires large 
gauge groups and small scalar representations. 


For such representations the potential in lowest- 
order perturbation theory, being restricted to re- 
normalizable couplings, is not sufficiently com- 
plicated to allow for complete spontaneous sym- 
metry breaking. 

We clearly have not exhausted all gauge groups 
and all scalar representations. One interesting 
possibility that we have not fully explored is that 
the symmetry breaking occurs beyond the tree ap- 
proximation, following Coleman and Weinberg.” 
As we noted above a perturbation-theory calcula- 
tion of the potential is unreliable for small values 
of the classical fields (which corresponds to small 
momenta), However, perturbationtheory is reliable 
for asymptotically free theories, for large values 
of the classical fields. Therefore there exists the 
possibility that the potential has a stable minimum 
for large vacuum expectation values of the scalar 
fields and that these are calculable using the re- 
normalization group. This possibility is rendered 
more likely by the fact that in many models the 
quartic scalar coupling constants (A and 7) can be 
negative for p? =—” but become positive as —p? in- 
creases, approaching zero from above as p? — — œ 
(see Fig. 13). Another possibility is to include 
Yukawa-like couplings between the scalars and the 
fermions. In the absence of gauge fields Yukawa 
couplings are never asymptotically free,?? but it is 
easy to see that this is no longer true in gauge the- 
ories. We have not explored this possibility in 
detail. 

In conclusion it appears to be very difficult, if 
not impossible, to construct theories which are 
both asymptotically free and which contain no 
massless vector mesons in perturbation theory. 
Such a model would be of great interest in provid- 
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NY, 


FIG. 13. The coupling-constant trajectories (a and 6 
move along one of the directed lines as ¢ increases). 
The gauge group is SU(6)xSU(6) and the scalars trans- 


form according to the (6x 6) representation. 


ing a theory which possesses an S matrix in per- 
turbation theory and is asymptotically free. How- 
ever, we do not believe that it would be of physical 
interest. First, if such a model exists it must be 
rather complicated and almost unique. Second, as 
we have seen above, in order to achieve asymptotic 
freedom it is necessary to incorporate many fer- 
mions and render f, very small. This might im- 
ply that 8, has an IR-stable fixed point at a small 
value of g*, thus restricting the physical coupling 
constant g(0, g) to be small if we are to remain 
within the domain of attraction of the origin. Such 
a model would probably not be useful as a theory 

of the strong interactions. Furthermore if p, is 
small, for small g, then the approach to scaling, 
according to Eq. (4.16), is very slow. This is 
clearly not a desirable feature. Therefore, to con- 
struct a physically meaningful model of the strong 
interactions one probably must pin one’s hopes on 
the possibility of dynamical symmetry breaking. 


VII. CONCLUSIONS 


The theories proposed in this paper are incom- 
plete. The main problem which requires investiga- 
tion is whether one can obtain an infrared sensible 
theory without explicit Higgs mesons. One might 
expect, on physical grounds, that the infrared sin- 
gularities induced by the gauge charges (color) are 
so strong that they must be completely shielded, 
so that only objects neutral under the gauge group 
could exist. This is an exciting possibility which 
might provide a mechanism for having a theory of 
quarks without real quark states. Whether this 
can be realized or whether the theory will exhibit 


dynamical symmetry breaking deserves much at- 
tention. 

What we have achieved so far is to find a large 
class of asymptotically free theories. We have 
shown that all semisimple gauge theories are in 
this class, as well as many theories involving fer- 
mions. We have explored the consequences of this 
asymptotic freedom with respect to deep-inelastic 
scattering and we have constructed some models 
which contain scalar mesons. Finally let us recall 
that the proposed theories appear to be uniquely 
singled out by nature, if one takes both the SLAC 
results and the renormalization-group approach 
to quantum field theory at face value. 


APPENDIX A: THE GAUGE DEPENDENCE OF THE 
RENORMALIZATION-GROUP EQUATIONS 


Let us consider the time-ordered product of 
gauge-invariant operators O, which are multipli- 
catively renormalizable. These could be, say, 
cauge~invariant currents. This renormalized 
Green’s function, which we denote by G'”, will sat- 


> 


isfy the renormalization-group equation: 


ð 9 5 
ue alg, eae v;(g, a)(1- 203°.) |e )=0, 
(A1) 


where we have used Eq. (2.21) and the sum runs 
over the anomalous dimensions of the operators 
O;. 
The derivative with respect to the physical gauge 
parameter can be eliminated using gauge invari- 
ance. This implies that the unrenormalized ampli- 
tude Gi”) is independent of the bare gauge parame- 
ter a,: 


ð 
Oa 


=, (A2) 


& y (A fixed 


Ce A) 


u 


When we recall Eq. (2.14) and express G{P in 
terms of the renormalized Green’s function we 
have 


E -ao(g, a) =+ elg, a) + Dols a) 


xG s, œ, u)=0, (A3) 


where 
Za 
= 4 
olg, a) Oy TE TN Se 
dln Z; 
o,(g, a)=Z, Jq Ae ’ (A5) 
u Sy: 1X 
Og 
(ga)=Z52| (A6) 
u Eu’ 1xe 


This equation expresses the fact that a change in 
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gauge parameter for gauge-invariant Green’s func- 
tions can be reabsorbed by a change in coupling 
constant and in the scale of the operators. It can 
be used to eliminate the 9/aqa in Eq. (A1). 

Finally let us note that to lowest order 6(g, a) is 
independent of a. This is essentially because, to 
lowest order, the coupling constant is unique. In 
other words if we change œ we might change g to 
g' however, g’ (being the value of the three-point 
function at some point) can be expressed as a pow- 
er series in g: 


g'=g+0(g°). (A7) 


The term of order g° might depend on a; however, 
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The Bjorken-Johnson-Low limit is used in a canonical current-algebra framework to estimate the 
magnitude of weak radiative corrections in second and fourth order of the weak coupling constants for 
unified gauge theories with strong interactions accounted nonperturbatively. We study, in particular, the 
question of violation of selection rules such as those for parity and strangeness ([AS| = 1 to second 
order and |AS| = 2 to fourth order) for a large class of models, in which strong interactions are 
generated by Abelian gauges (the neutral-vector-gluon theory) or a special set of non-Abelian gauges. 


1. INTRODUCTION 


Unified gauge theories of weak and electromag- 
netic interactions are threatened by the possibility 
that they may, in general, lead to intolerable mag- 
nitudes for the violation of selection rules such as 
those for parity and strangeness via radiative 
corrections. The root of the problem rests with 
that appealing feature of the framework which 
renders possible the description of a world of 
interactions having radically different strengths 
at low energies (momentum transfers) in terms 
of the same coupling constants (g, g’, ...). Within 
such a framework a hierarchy of “effective 
strengths” emerges for the electromagnetic and 
the weak interactions at low energies in the tree 


approximation (i.e., a=g°/47 for the electromag- 
netic and Gm? =g°m?/u? for the weak; m is some 
typical hadronic mass and uy is the typical mass 

of a heavy weak gauge boson) as a consequence of 
the large difference of mass between the photon 
and the weak intermediate vector mesons. How- 
ever, this hierarchy, which also ceases to exist 
at high energies, is not generally preserved by the 
weak radiative corrections. The obvious reason 

is that loop integrations appear which are either 
divergent or slowly convergent and for which there 
is not a significant numerical suppression cor- 
responding to the large intermediate boson mass. 
Such a suppression would indeed only be present 
where loop integrations converge fast enough to 
allow one to neglect, to a good approximation, the 


